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Abstract 

The dramatic size evolution of early-type galaxies from z ~ 2 to poses a new challenge in the theory 
of galaxy formation, which may not be explained by the standard picture. It is shown here that the size 
evolution can be explained if the non-baryonic cold dark matter is composed of compact objects having 
a mass scale of ~ lO^M©. This form of dark matter is consistent with or only weakly constrained by the 
currently available observations. The kinetic energy of the dark compact objects is transferred to stars 
by dynamical friction, and stars around the effective radius are pushed out to larger radii, resulting in a 
pure size evolution. This scenario has several good properties to explain the observations, including the 
ubiquitous nature of size evolution and faster disappearance of higher density galaxies. 
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1. Introduction 

It is well established that early-type galaxies have 
mostly been built up before redshift z ^0.8 (e.g., Cirasuolo 
et al. 2007; Pozzeti et al. 2007). However, recent obser- 
vations have revealed that passively evolving, early-type 
galaxies at higher redshifts (z~ 1~2) are surprisingly much 
more compact than the likely present-day descendants of 
the same stellar mass (e.g., Longhetti et al. 2007; Trujillo 
et al. 2007; Cimatti et al. 2008; Damjanov et al. 2009, 
hereafter D09). Such compact galaxies are not found in 
the local universe, and they must disappear (decrease in 
comoving number density by a factor >5000, Taylor et al. 
2009) by structural evolution, meaning that the evolution 
must occur ubiquitously, rather than stochastically, for all 
such galaxies. In the plane of M^:-re (total stellar mass and 
the effective radius including half of the projected light), 
they must be "puffed up" with a typical size increase by a 
factor of 2-3 (a factor of more than 10 in stellar density) to 
move into the region populated by the present-day galax- 
ies, without significant increase in stellar mass (Fig. 1). 
Furthermore, this evolutionary process seems to work in 
a broad context of galaxy formation spanning wide ranges 
of mass, radius, and redshifts (D09). 

This size evolution problem thus presents a new major 
challenge for the theory of galaxy formation and evolution. 
Several scenarios have been proposed, but none of them 
has convincingly explained it. Here, a new scenario is pro- 
posed that may contribute to solve the problem. It is as- 
sumed that the dominant component of the non-baryonic 
cold dark matter in the universe is dark compact objects 
(DCOs), like primordial black holes (PBHs) or substruc- 
tures in dark haloes, having a mass scale of M, 10'^ M©. 
This possibility is not yet excluded by currently available 
observations (sec §4). 

In §2, previously proposed scenarios for this problem 



and their potential difficulties are briefly reviewed. The 
new scenario of this work will be described in §3, followed 
by some discussions in §4 and §5. 

2. Proposed Scenarios for the Size Evolution 

Dissipationless "dry" mergers of old stellar systems can 
increase galaxy sizes without new star formation activ- 
ity, but stellar mass should also be increased. Bezanson 
et al. (2009) argued that mass increase with a scaling 
of A/* cx Te, which is expected in the case of major or 
equal-mass mergers, cannot solve the problem because it 
would violate the constraint from the stellar mass found 
in present-day early- type galaxies. 

Minor dry mergers (re (x in the limit of very large 
mass ratio) may be more helpful to solve the problem, 
as demonstrated by a simulation of Naab et al. (2009). 
However, their simulation is only for one galaxy, and it 
must be further studied whether this scenario can explain 
the size evolution in the general context of galaxy for- 
mation (D09; Williams et al. 2009). All high-z compact 
galaxies must experience many minor mergers to gain sub- 
stantial mass and size increase, and they would have to 
be all dry to keep the old stellar population, which seems 
somewhat contrived. It should be noted that the largest 
mass (~ IO^^Mq) of z ^ 1 galaxies in the D09 sample is 
already comparable with that of the local SDSS sample 
(see Fig. 1), and further mass increase may contradict 
the data even in the case of minor mergers. Theoretically, 
minor dry mergers should be effective for massive and 
high redshift galaxies (Naab et al. 2009; Hopkins et al. 
2009a) , but the D09 data indicate that the required degree 
of size evolution is similar for galaxies in the mass range of 
1O^"'^-1O^^M0, and lower mass galaxies disappear 
faster at higher redshifts (Fig. 1). Nipoti et al. (2009) 
have argued that constraints from gravitational lensing 
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Fig. 1. The effective radius (re) versus total stellar mass (M«) and stellar density (p*e, the moan within r^) of early-type galaxies. 
The red (filled) and blue (open) symbols are for high-rcdshift galaxies (1.46 < z <2 and 1.1 < z < 1.46, respectively), derived or 
complied by D09. Small grey dots are mass estimates by D09 for the local SDSS early-type galaxy sample. All the stellar mass 
data are corrected to a common stellar initial mass function (IMF) by D09. The contours of the energy transfer time scale a^re 
depicted assuming M, = 10^ Mq for DCOs in the both panels (ti/2 ocM^^). The arrows indicate the evolutionary directions expected 
by a pure-size evolution (magenta), ideal minor mergers (green, oc M^), and major (equal-mass) mergers (black, oc M«). 



do not allow size evolution by more than a factor ^^1.8 by 
dry mergers, which seems too small to completely solve 
the size evolution problem. 

Puff up of galaxies by mass loss is another possibility, 
such as stellar mass loss (D09) or by quasar feedback (Fan 
et al. 2008). However, stellar mass loss is not satisfactory 
because of the small amount of mass loss and short time 
scales (D09). The physical process of quasar feedback is 
rather uncertain, and it seems difficult to puff up quiescent 
galaxies in passive evolution phase after they lost most of 
interstellar gas (Bezanson et al. 2009). It may require a 
significant fine-tuning to explain the ubiquitous size evo- 
lution by rather stochastic quasar feedback processes. 

3. Size Evolution by DCOs 

Various observations of nearby early-type galaxies 
based on stellar kinetics (e.g., Cappellari et al. 2006), X- 
ray emitting interstellar gas (e.g., Fukazawa et al. 2006), 
and strong/weak gravitational lensing (e.g., Gavazzi et 
al. 2007) have shown that stellar mass is dominant at 
inner radii while dark matter dominates at outer radii, 
and the transition occurs at around the effective radius. 
Therefore, if there is significant transfer of kinetic energy 
from dark matter to stars, it would puff up stars at r ^ re 
resulting in a pure-size evolution of the stellar component. 

This energy transfer occurs by dynamical friction. 
When stars have a stellar mass density p^e and the 
Maxwellian velocity distribution with one-dimensional 
dispersion cr,e, the dynamical friction timcscale for a DCO 
having a velocity V is given as (Binney & Tremaine 1987): 



V 



V 47rG21nAp*eM.B(X) 



(1) 



where the function B{X) = erf(X) -2Xe~^' and 
X = V/{V2a^e)- The Coulomb logarithm In A = 
ln[6max^^/(GM.)] becomes 13.74 and 7.60 if we take the 
maximum impact parameter fomax to be 10 kpc (the galaxy 
scale) and {p^e/M,)^^/^ (the mean separation of DCOs), 
respectively, for typical values of M, = IO^Mq, V = 200 
km/s, and pn-e = 10^°Mokpc~^. In the following InA = 10 
is adopted. 

Since ^df is strongly dependent on V, we should care- 
fully take the effective mean about V. The total kinetic 
energy i;cx / {M,V'^ /2)f{V)d^V of DCOs wiU decrease by 
a factor of 1/2 in a time of ti/2 = E/ {2E), where f{V) is 
the phase space distribution function of DCOs. Assuming 
that the initial velocity distribution of DCOs is the same 
as that of stars, ti/2 is found as: 



tl/2 = 



/372)r(5/2) 



3V6/o°°B(V^)e-^da; 



O.BHa 



(2) 



where tav is t^i when V is the mean value (T/^)^/^ — \/3a^e. 

To compare with the observed data in Fig. 1, we con- 
vert CT»e into stellar mass using the direct correlation be- 
tween the two (i.e., equivalent to the Faber- Jackson re- 
lation for a constant stellar mass-to- light ratio). The 
virial relation a^^ ~ GM^/r^ is another option of this 
conversion, but we take the former because the M*-cr*e 
correlation is tight and the latter could be affected by 
the size evolution (as discussed here) and contribution 
from dark matter. We expect rather mild evolution of 
gainst the size evolution, considering the roughly 
flat velocity profiles of elliptical galaxies including dark 
matter (e.g., Mamon & Lokas 2005). Using observed ve- 
locity dispersions^ CTobs of the SDSS early-type galaxies 



The Bernardi et al.'s (Tobs has been corrected to an aperture 
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(Bernardi et al. 2003) and the stellar mass estimates of 
D09, A/* =4.7-lO"(cr*e/2OOkms-i)3-98M0 is found, and 
hence the energy transfer time scale becomes 



A/* 


-0.247 


re 
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M, 






1 kpc 




105A/q_ 



(3) 



Here, the stellar mass fraction of 0.42 within (corre- 
sponding to the Scrsic index n = 4, e.g., Mamon & Lokas 
2005) has been used to convert p*e (the mean density 
within Tg) into and r^. 

This time scale is compared with the observed data of 
early-type galaxies in the plane of M^-r^ as well as p*e-fe 
(Fig. 1). The high density galaxies (p*e ^ lO^^Afokpc"^) 
of the D09 sample is found only in the higher rcdshift 
bin oi z > 1.46 (red symbols in Fig. 1), indicating that 
these galaxies should disappear on a shorter time scale 
of a few Gyr, while galaxies of lower densities (p*e ^ 
lO^''Af0kpc~'^, blue symbols) seem to disappear more 
slowly on a time scale of ~ 10 Gyr. Fig. 10 of DG9 has 
shown this trend even more clearly in a p*e--z plot, where 
galaxies having higher stellar density disappear at higher 
redshifts. This trend is expected in the proposed scenario, 
because galaxies having higher stellar density are located 
in the region of shorter ti/2, as indicated in Fig. 1. 

The ii/2 of high-z early- type galaxies with the largest 
mass is larger than the age of the universe. We may need a 
larger DCO mass than the assumed 10^ Mq, or some other 
mechanisms, for these galaxies to evolve into the location 
of the local galaxies. On the other hand, Mancini et al. 
(2009) has recently reported that the sizes of the most 
massive early- type galaxies at z ^ 1.5 are similar to those 
of local galaxies, indicating a possibility that the size shift 
in the Mt.-re plane is smaller for more massive galaxies. 
In this case DCOs with A/, ^ IO^A/q could be sufficient 
to solve the problem. 

It should be noted that stellar density is larger than that 
of dark matter at r ^ Tg (e.g., Mamon & Lokas 2005), and 
we do not expect a strong structural evolution of stars 
at r ^ Tg.^ This is in contrast to the change of stellar 
distribution by stellar remnant black holes that have the 
same initial density profile as that of stars (e.g., Merritt 
et al. 2004). The small effect at r ^ is consistent 
with the observations indicating that the core densities 
of high-z elliptical galaxies are not much different from 
those of local galaxies (Bezanson et al. 2009; Hopkins et 
al. 2009b). The expected mild evolution of velocity dis- 
persion is also consistent with the observational result of 
Cenarro & Trujillo (2009).^ 

radius of re/8, and the central velocity dispersion is typically 
~ 1.4 times higher than that around (e.g., Cappellari et al. 
2006). Therefore o"obs = 1.4cr«e has been adopted here. 
^ The accumulation of DCOs in a galactic center may affect the 
innermost stellar distribution, such as stellar core formation, but 
it would also be complicated by the existence of super-massive 
black holes. This is an interesting issue to be investigated, but 
beyond the scope of this letter. 

But see also van Dokkum et al. (2009) who have reported a 
large velocity dispersion for one galaxy at z = 2.186. Velocity 
dispersion measurements of early- type galaxies at 2 ^ 1 arc still a 
difficult task, and more observational studies are highly desired. 



To conclude, the proposed scenario has several features 
that are good to explain the observed trends: (1) the evo- 
lution should be pure-size evolution, (2) it works ubiqui- 
tously for all galaxies in a non-stochastic manner on cos- 
mological time scales, (3) higher density galaxies should 
disappear on shorter time scales, (4) stellar density at 
r <^re is not much affected, and (5) evolution of velocity 
dispersion should be mild. 

4. Observational Constraints on DCOs 

There arc no strong observational constraints to exclude 
the possibility of lO^Af© DCOs as the dominant form of 
dark matter (for reviews see Mack et al. 2007; Ricotti et 
al. 2008, hereafter M07 and R08). DCOs of < lO^Af© 
are marginally allowed by constraints from strong lensing 
of compact radio galaxies (Wilkinson et al. 2001) and 
millilcnsing of gamma-ray bursts (Nemiroff et al. 2001). 
On the other hand, anomalous flux ratios in gravitational 
lenses may be explained by DCOs of lO^~^A/0 rather than 
dark matter of particle-mass scale (Mao et al. 2004). 

Yoo et al. (2004) declared "the end of MACHO era" 
by closing the window of Af, ~ 30-10^ A'/q based on the 
existence of wide binary stars in our Galaxy halo, but 
the robustness of their conclusion has been questioned by 
Quinn et al. (2009) using new data. 

DCOs of Af, < IQi^Mq do not affect the stellar distribu- 
tion of dwarf spheroidal galaxies (dSphs) like Draco in the 
Milky Way halo (Jin et al. 2005), and lO^-'^A^© DCOs 
may even be helpful to solve some problems in galaxy for- 
mation such as disk heating (Lacey & Ostriker 1985) and 
the core/cusp problem (Jin et al. 2005). Sanchez-Salcedo 
& Lora (2007) has claimed that IQ^Mq DCOs are ex- 
cluded as dark matter from the survival of the kinetically 
cold stellar clump found in the Ursa Minor (Klcyna et al. 
2003, hereafter K03). However, quantitative constraint on 
the velocity dispersion of the cold clump was not derived 
by K03, making a quantitative analysis difficult. The ori- 
gin and history of the clump is also still uncertain (see, 
e.g., Mufioz et al. 2008; Penarrubiaet al. 2009). It should 
be noted that UMi is showing substantial morphological 
distortions, unlike most other dSphs (K03). 

Afshordi et al. (2003) excluded DCOs heavier than a 
few times lO^Af© as dark matter, using the small-scale 
density ffuctuations inferred by Lya forests. However, 
this analysis is depending on a modeling of complicated 
baryonic physics about Lya forests that may include sig- 
nificant systematic uncertainties. 

Mil & Totani (2005) and MapeUi et al. (2006) discussed 
the expected number of X-ray sources in nearby galaxies 
considering radiatively inefficient accretion flow onto halo 
black holes of A/, ^ lO'^A/©. R08 set strong constraints on 
PBHs as dark matter by requiring that the accretion ac- 
tivity onto PBHs does not affect the anisotropy of the cos- 
mic microwave background. However, these calculations 
include many assumptions and large uncertainties espe- 
cially about accretion physics and feedback effects. R08 
discussed about angular momentum estimating it from the 
normal cosmic density fluctuation. However, the proper 
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mean separation of 10^ A/0 dark matter DCOs at z 1000 
is about 10 pc, which is comparable to the Bondi radius 
at that epoch estimated by R08. Then interaction with 
nearby DCOs would result in significant angular momen- 
tum that should prevent efheient accretion. 

5. Discussion 

Although DCOs are less popular than elementary par- 
ticles as the candidate of dark matter, several formation 
scenarios have been proposed (see M07 and R08 for re- 
views). PBHs can be formed by density fluctuation of a 
modest amplitude or softening of equation of state, when 
the fluctuation scale enters the event horizon of the uni- 
verse. The mass scale of IO^Mq corresponds to the hori- 
zon mass (Mhor = tc^/G) at a temperature of ~ MeV, i.e., 
close to the annihilation era. Note that only a small 
portion 10"^) of the universe should be converted into 
PBHs to explain the dark matter density, and hence it 
may not necessarily affect the primordial nucleosynthesis. 

DCOs need not to be very compact like PBHs, but sub- 
structures in dark haloes generally found in simulations of 
structure formation should also contribute to dynamical 
heating of stars. However, the mass fraction of substruc- 
ture must be of order unity to explain the size evolution, 
which seems rather unlikely within the standard picture 
because less than ~10% of dark mass within the virial ra- 
dius (and even smaller fraction at inner radii) is found in 
substructures in recent simulations (Springel et al. 2008). 

The black hole masses found in active galactic nuclei are 
always ^ IO^Mq. Black holes heavier than IO^Mq have 
been discovered in quasars beyond z > 6 (e.g., Willott et 
al. 2003), and their existence at such high redshifts is not 
easy to explain by evolution starting from stellar mass 
black holes (e.g., Tanaka & Haiman 2009). It is obvious 
that the explanation would become easier, if dark matter 
is composed of lO^A/0 DCOs. 

More realistic, numerical studies of galaxy size evo- 
lution would be required to conflrm this hypothesis. 
Observational searches for DCOs by various approaches 
(as discussed in §4) may detect DCOs in the future. 
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"The Next Generation of Physics, Spun from Universality 
and Emergence" from the Ministry of Education, Culture, 
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